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Where to find a dark matter sterile neutrino? 
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We propose a strategy of how to look for dark matter (DM) particles possessing a radiative 
decay channel and derive constraints on their parameters from observations of X-rays from our own 
Galaxy and its dwarf satellites. When applied to the sterile neutrinos in keV mass range, it allows 
a significant improvement of restrictions to its parameters, as compared with previous works. 
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Introduction. It was noticed long ago that a ster- 
ile neutrino with the mass in the keV range appears to 
be a viable DM candidate |jyj. Moreover, being a warm 
DM, sterile neutrino ease the tension between observa- 
tions and predictions of the cold DM model on small 
scales. The interest to this scenario is revitalized since 
the discovery of neutrino oscillations (see e.g. for re- 
view). Indeed, one of the simplest ways to explain this 
data is to add to the Standard Model several (at least 
two) gauge singlet fermions - right handed, or sterile 
neutrinos. It has been demonstrated recently Q that 
a simple extension of the Standard Model by three sin- 
glet fermions with masses smaller than the electroweak 
scale, dubbed the i/MSM in 3], allows to accommodate 
the data on neutrino masses and mixings (with the ex- 
ception of the LSND anomaly), provides a candidate for 
DM particle in the form of sterile neutrino, and allows to 
explain baryon asymmetry of the Universe. The simplic- 
ity of the model, the similarity of its quark and lepton 
right-handed sectors, together with a considerable num- 
ber of phenomena it can simultaneously describe, forces 
us to take this model seriously and thus provides addi- 
tional motivations for study of keV mass range sterile 
neutrino as a DM candidate. 

The sterile neutrino has a radiative decay channel, 
emitting a photon with energy E — ms/2 {nis being the 
mass of sterile neutrino) . Parametrically, the decay width 
is proportional to sin^ 29 Q , where 6 is the mixing 
angle between active and sterile neutrino. 

If such a neutrino is a main ingredient of the DM, it is 
potentially detectable in various X-ray observations. The 
most obvious candidates include diffuse extragalactic X- 
ray background (XRB) 5, 6, 7,'8'|; clusters of galaxies IB, 
0, ilfl] ; galaxies ^\ , including our own. 

The aim of the present Letter is to discuss the best 
strategy to search for a DM sterile neutrino and to derive 
the constraints on its properties. Although we concen- 
trate on sterile neutrino, the constraints we get can be 
applied to any DM candidate with a radiative two-body 
decay channel in a keV range. We analyze various types 
of astrophysical objects and show that the strongest con- 
straints on sterile neutrino are coming from neutrino de- 
cays in the Milky Way halo and in particular in the halo 
dwarf galaxies. These objects were not considered previ- 
ously in this context. The existing XMM Newton and 



HEAO-1 data allow us to improve over previous con- 
straints, highlighting the potential of new optimal cites 
for the searches of the signal from the sterile neutrino 
decay. 

Dark matter halo of the Milky Way (MW). En- 
ergy flux produced by the DM decay from a given direc- 
tion into a solid angle Qfov <SC 1 is given by 
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where F is the radiative decay width of sterile neutrino. 
To determine the MW contribution into the flux one 
needs to know the distribution of the DM in the halo. 

Mass distribution within MW has been modeled by 
many authors. Characteristics which are relevant for 
our study are tightly constrained by the wealth of de- 
tailed data available for this system. As a reference we 
will choose the mass distribution derived recently in , 
where physically interesting models were selected by im- 
posing additional constraints based on a theory of halo 
formation. At large r the halo density can be described 
by the Navarro- Frenk- White (NFW) profile /3nfw('') — 
¥^—! — The MW halo parameters of favorite mod- 
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els obtained in correspond to Afvir — 1-0 x lO^^M©, 
= 21.5 kpc and numerical constant a ~ 1.64. 
In the region of r relevant for our study the halo density 

can be also approximated by the isothermal profile 
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where Vh corresponds to contribution of the DM halo 
into the Galactic rotation curve in its flat part, w/j « 170 
km/s, see e.g. [llj |. 

The NFW density profiles and (O produce identical 
fluxes from the direction of the Galactic anti-center if 
Tc w 4 kpc, giving J^phaiodr w 0.7 x 10^2 GeV/cm^. 
They also follow closely each other in the range 3 kpc < 
r < 80 kpc (difference being less than 5%). Therefore, in 
estimates of the flux from directions which are outside of 
20° circle around the Galactic center both NFW and 
give the same results. 

The halo density profile is less certain in the region 
r < tq. In two distinct types of models were con- 
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sidered, with and without exchange of angular momen- 
tum between DM and baryons. In the model without 
momentum exchange the DM density profile at r < 10 
kpc diverges even faster than NFW profile. In the model 
with angular momentum exchange the DM density pro- 
file at 2 kpc < r < 10 kpc is less singular and rather 
resembles the isothermal sphere Eq. ^ with rc « 4 kpc, 
but the DM density is larger at r < 2 kpc as compared 
to the isothermal sphere. Therefore, one can use Eq. ^ 
with Tc = 4 kpc as a lower limit on the DM density 
when calculating fluxes from DM decays, and, therefore, 
for putting a conservative bound on the sterile neutrino 
mixing angle (e.g. the halo density at the Sun position 
in the model Eq. ^ is 0.25 GcV/cm , which is smaller 
than the accepted value 0.3 GeV/cm'^ E3)- Utilizing the 
NFW profile at all r can only strengthen the bounds. 

In the model Eq. the DM fiux from the direction 
(6, 1) (in galactic coordinates) into the solid angle fifov ^ 
1, measured by an observer on Earth is given by 
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and cos d — 



where Lo ^ 

cos 6 cosZ. For example, -Fmw(90°)/Fmw(180°) ~ 1.52. 

Search for a preferred observation. Let us com- 
pare the Galaxy contribution to the DM decay flux, com- 
puted above, with those of other astrophysical objects. 

(i) XRB. Although the DM has a very narrow radiative 
decay width, the cosmological DM decay contribution to 
the XRB gets significantly broaden due to the contribu- 
tions from various red shifts. The resulting differential 
flux for E < TOs/2 is given by 
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where ~ 1.2 x lO^^GeV/cm'^ is the average DM 
density in the Universe, J^a and flm are the cosmological 
constant and matter contributions to the density of the 
Universe. The paper 5] looked at restrictions from XRB, 
assuming that DM is uniform up to very small z. This 
question was further addressed in Q and finally the 
most stringent constraint in the (sin^ 20, Wg) plane from 
XRB was obtained recently in |^. 

One can compare flux @, integrated over all E with 
the galactic contribution: Fmw/^xrb = T^TcHq, where 
~ Pmw/Pdm ^ l*-*^ overdensity of the Galaxy as 
compared to the average density of DM in the universe. 
With Tc ^ 4 kpc, one arrives to the conclusion that 
^Mw/^xRB 1, i-e. Galactic contribution is compara- 
ble with the total DM decay flux from all red shifts. 

However, for a modern X-ray instrument with good 
spectral resolution AE ^ nis (e.g. XMM-Newton) one 
should compare contributions from the Galaxy and from 
uniform cosmological distribution into XRB within AE. 
The ratio Fmw/Fxrb then gets enhanced by the factor 



E/ AE which is 10 50 for EPIC cameras on board of 
XMM-Newton, i.e. for XMM the Galactic signal is 1 to 2 
orders of magnitude stronger than the contribution from 
uniform distribution of DM in the universe. 

(ii) Clusters. Let us now analyze the flux enhancement 
from Coma and Virgo clusters of galaxies Ii|l3. The 
DM distribution in these clusters can be approximated 
by isothermal 13-model 13, 14] 
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Integrating it according to Eq. |^ and using the fact 
that the Coma cluster is located perpendicularly to the 
Galactic plane, and has core radius rdustcr — 0.3 Mpc 
and central density Pcoma = GeV/cm'^ (see e.g. [l^ 
for discussion), we get -FMw(90°)/Fcoma ~ 0.25. Simi- 
lar estimate for the center of the Virgo cluster (pvirgo — 
2.3 GeV/cm'^, rdustcr ~ 10 kpc) shows that Galactic con- 
tribution is ~ 10%. Therefore, one could conclude that 
clusters are preferable objects for DM detection 6,, A]. 

These results have been recently reanalyzed by [lOl| . 
where it was shown that the actual restriction from clus- 
ters is only by about a factor of 3-4 (in 9'^ for given rris) 
better than those from XRB. The reason for such a mod- 
est improvement is the following. The detection of the 
DM decay line in galaxy clusters is complicated by the 
fact that most of them show strong emission precisely 
in keV range. Indeed, the virial theorem immediately 
tells us that the temperature of the intercluster medium 
is Tgas ~ Gn rUpTZp'^f^id'^ , where d is characteristic size. 
For overdensity TZ 10'^ and size d ~ 1 Mpc, the tem- 
perature Tgas is always in the keV range, which makes it 
hard to detect a DM decay line against a strong X-ray 
continuum. 

(in) Dwarf galaxies. There should be an enhance- 
ment of the flux in directions of dwarf spheroidal galax- 
ies, which are satellites of the MW. Promising satellites 
with large mass to light ratio are Draco and Ursa Minor. 
Density proflles of both galaxies can be modeled by the 
isothermal sphere with Vh ~ 22 km/s and Tc ~ 100 pc 
[Tsll . This gives for the contribution to the flux from the 
dwarf galaxy (along the line which passes through the 
core of the sateflite) /pdr « 3.3 x 10^^ GeV/cm^. 

Contribution of the Galaxy halo flux in directions of 
both sateUites « 1.0 x 10^^ GeV/cm^. Therefore, in di- 
rections of both satellites the 4-fold local enhancement 
of the flux is expected, while the total flux matches the 
flux from the central region of the MW and the flux from 
clusters. The advantage of observing the dwarf satellites, 
as compared to clusters or to the Galactic center, is lower 
level of X-ray background contamination and clear signa- 
ture of the signal, namely, local flux enhancement within 
single fleld of view of X-ray telescope. 

Therefore, we see that the local DM halo (especially 
dwarf satellite galaxies) can provide the strongest restric- 
tion on the parameters of the sterile neutrino as DM can- 
didates. 

Restrictions from local halo DM contribution. 
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FIG. 1: Upper panel: Method used to obtain restrictions on 
DM decay in MW from blank sky XRB observations. The 
data is fitted by a power law (reduced ~ 1-07 for 82 d.o.f.) 
and XSPEC vll.3.2 is used to put a 3a limit on the presence 
of DM line (via command "error <line norm> 9.0"). 
Lower Panel: Method used to obtain restrictions from LMC. 
Flux rapidly decreases for E > 2 keV, most of the data points 
at higher energies are zero within statistical uncertainty. The 
solid green line is the sum of the total flux plus 3a per energy 
bin. At i5 > 2 keV the bound is dominated by errors and 
therefore can be improved significantly by increasing statis- 
tics. 
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FIG. 2: 3a restrictions on DM decay line from Milky way halo 
(XMM and HEAO-1 observations). Dashed line: total flux + 
3a restriction from LMC. Also shown the previous strongest 
limit from the clusters of galaxies [lo|. 

We have analyzed XMM blank sky observations of [l^ 
(exposure time ~ 200 ksec), taking into account MW 
contribution and putting a 3cr bound on the DM line 
flux in the region 1 keV < E < 7 keV (FIG. UJ. The 
exclusion plot in the region nis > 6 keV was obtained 
by using the HEAO-1 measurements of XRB 0|. As 
spectral resolution of HEAO-1 is about 25%, the resulting 
correction due to MW is not so drastic as compared to 



the results of Ref. . 

Unfortunately, X-ray observations of Draco and Ursa 
Minor dwarfs are not available currently. However, ap- 
proximately the same signal is expected from the Large 
Magellanic Clouds (LMC) albeit with larger uncertain- 
ties. Using isothermal sphere, with the following halo 
parameters, Vh ~ 50 km/s and Tc ~ 1 kpc p^ . we obtain 
J pdr Ri 2.8 X 10^2 GeV/cm^ for the LMC contribution, 
while MW halo contribution in the direction of LMC is 
again w 1.0 x 10^^ GeV/cm^ (the use of the NEW profile 
gives even larger value for flux). 

Therefore, as an illustration we have processed one of 
the observations of the LMC (XMM obs ID: 0127720201, 
exposure ^ 20 ksec). Subtracting the blank sky back- 
ground one sees that the flux is zero within statis- 
tical uncertainty for E >2 keV. Similar reduction of the 
background is expected for dwarf satellites. In this situ- 
ation we put an upper limit on the flux of the DM decay 
line, by demanding it to be smaller than total flux plus 
its 3f7 error in an energy bin equal to spectral resolu- 
tion, see Fig^ One can see that increasing the exposure 
time, one can signiflcantly lower the restriction especially 
in the region E >2 keV. 

Discussion. In this paper, we have shown that the 
best objects for the search of the DM with radiative de- 
cay channel is the MW halo, including dwarf satellite 
galaxies (e.g. Draco and Ursa Minor). To illustrate this, 
we put the strongest restrictions on parameters of sterile 
neutrino (i) searching for the MW DM decay signal in 
the blank sky XRB; (ii) using LMC observations. One 
can see that improving statistics one can significantly im- 
prove bounds from such objects. Our analysis also shows 
that improvement of the spectral resolution of X-ray in- 
struments (even by means of decreasing imaging capa- 
bilities) is crucial to continue the search of DM decay 
line. 

Of course, all constraints on the sterile neutrino mixing 
angle, derived from X-ray observations, suffer from the 
uncertainties in the DM profiles. However, in our analy- 
sis we tried to be as conservative as possible. Namely, we 
present only the results coming from the study of direc- 
tions which are away of the Galactic center. In this case 
the most relevant parameter, which infiuences the line of 
sight integral Eq. is Vh, corresponding to the contri- 
bution of the DM halo to the Galactic rotation curve at 
large distances, away form the core. The total (DM plus 
baryons) rotational curve is measured directly. Although 
the subtraction of the baryon contribution is model de- 
pendent, this dependence is weakest for the directions 
opposite to the Galactic center. 

Moreover, we have found, quite remarkably, that the 
line of sight integral Eq. ^ is roughly the same for all 
studied DM dominated objects, from cosmological back- 
ground to clusters of galaxies to dwarfs satellites. This 
suggests that the expected DM decay signal should be 
roughly the same for all of them. Better constrains 
are obtained for the objects whose X-ray background is 
lower. This makes dwarf satellites more suitable as com- 
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pared to clusters. In addition, in the most DM domi- 
nated satellites, such as Draco and Ursa Minor, not only 
the X-ray background is lower, but uncertainties due to 
subtraction of baryonic component from the galactic ro- 
tational curve are also minimized. The restrictions, based 
on the different (types of) objects, further minimize these 
uncertainties. 

The limits we derived here allow to strengthen the 
bounds on sterile neutrinos in different models of par- 
ticle physics, shed light on the possible mechanisms of 
their production in the early universe, and constrain dif- 
ferent astrophysical phenomena that might be related to 
their existence. 

(i) In the Standard Model with addition of just one 
sterile neutrino (assuming the absence of sterile neutri- 
nos above the temperature ^ 1 GeV and charge neutral- 
ity of the plasma) the relic abundance of sterile neutri- 
nos can be expressed through nis and 9 . This relation 
(quite uncertain, since the sterile neutrinos are mainly 
produced at temperatures 0(150) MeV, exactly where 
the description of the strongly interacting plasma is most 
complicated pol |') allows one to find, potentially, an up- 
per limit on the sterile neutrino mass in this particular 
scenario. Taking as a rough estimate the computation 
of and our X-ray constraints we arrive to an upper 
bound nis < 3 keV. This may be contrasted with the 
lower bound on the mass of sterile neutrino (derived in 
the same model with the same assumptions) coming from 
the analysis of the Ly-a forest data '2l'| : nis > 2 keV , 
TUs > 1.7 keV ,23j] leaving a very limited allowed mass 
range for a sterile neutrino. If a more recent result of [23 |. 
rUs > 14 keV is proven to be correct, and uncertainties re- 
lated to poor knowledge of QCD at the relevant range of 
temperatures happen to be not substantial, this scenario 



will be ruled out by cosmological and astrophysical obser- 
vations. The same conclusion is true ||25j for the z/MSM 
if the same assumptions about the initial state are taken 
for granted. This would make the production mecha- 
nisms of the sterile neutrinos related to inflation [2^ or 
large lepton asymmetries [2^ more important. 

(ii) The X-ray observations allow to predict the masses 
of active neutrinos in the z^MSM. In was shown, 
that the XRB limits of Q imply that the lightest ac- 
tive neutrino must have a mass < 3 x 10~^, provided 
nis > 1-8 keV. In case of normal hierarchy two other 
masses are given by the observed mass square differences 
y^A^^ 2, (8.5 - 9.5) X 10-3 eV and ^A^^ ^ 
0.04 — 0.06 eV. (In the case of inverted hierarchy masses 
of both neutrinos are \/ Amatm)- With the improved 
bound, derived in this paper, this result is true for even 
smaller sterile neutrino masses, rUs > 1.3 keV. 

(iii) Our constraints, combined with those of put 
severe bounds on explanation of pulsar kick velocities by 
dark matter sterile neutrinos, proposed in |,28i] and on 
the mechanism of early reionization by radiative decays 
of sterile neutrinos, suggested in |2^. 
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Note added. After finishing this paper we became 
aware of an independent analysis, where the possibility 
of measuring decaying DM particles from our Galaxy, 
using Chandra blank sky data, was considered [l^l . The 
conclusions, reached in that paper on MW, are similar to 
ours. 
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